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Non-unitary quantum mechanics has been used in the past to study irreversibility, dissipation 
and decay in a variety of physical systems. In this letter, we propose a general scheme to deal with 
systems governed by non-Hermitian Hamiltonians. We argue that the Schwinger-Keldysh formalism 
gives a natural description for those problems. To elucidate the method, we study a simple model 
inspired by mesoscopic physics -an asymmetric junction. The system is governed by a non-Hermitian 
Hamiltonian which captures essential aspects of irreversibility. 



Non-Hermitian formulations of quantum mechanical 
problems have attracted substantial interest in almost 
all areas of physics. In most cases, non-unitary ap- 
proaches are utiHzed to describe irreversible processes, 
-such as decay and dissipation-, in open quantum sys- 
tems. This type of problems have been addressed since 
the early days of quantum mechanics, when the com- 
plex eigenvalue method was pioneered to describe the 
a-decay of nuclei [J, Q|- Since then, non-unitary ap- 
proaches have been applied to theories of K and B me- 
son decay 0, |3|, scattering and absorption of particles 
by nuclei [J, nuclear reactions 0, 0| , multi-photon ion- 
ization of atoms 1 8, 9], optical resonators [10|| and free- 
electron lasers [ll|. Growth models have been investi- 
gated by mapping a master equation into a Shrodinger 
equation with a non-Hermitian Hamiltonian In 
the last decade, non-Hermitian theories have also been 
applied to condensed matter systems with localization- 
declocalization transitions. Pinning / depin ning of flux 
lines from columnar defects has been studied [l3| in type- 
II superconductors, and the effects of a single colum- 
nar defect on fluctuating flux lines has been considered 
by mapping to a non-Hermitian Luttinger Liquid model 
and usin g t he density matrix renormalization group 
(DMRG) GeneraHzed DMRG has been applied also 
to a one-dimensional reaction-diffusion model with non- 
unitary time evolution (l5| and a non-Hermitian spin- 
1/2 Heisenberg chain |16i]. The closing of the Mott 
gap has been investigated in the non-Hermitian Hub- 
bard model [l3| and evaporatively cooled Bose-Einstein 
condensates (BECs) have been studied by a non-unitary 
quantum dynamics approach . While all of the above 
studies used a variety of case-speciflc methods and were 
mainly concerned with either time-independent or tran- 
sient behaviors, in this letter we present a general frame- 
work, that we apply to a steady-state formulation of an 
irreversible system. 

To deflne our approach we start with the Schwinger- 
Keldysh (SK) formalism 3, 23|; which turns out to be 
very natural in the non-Hermitian case. In the SK for- 
mulation, the evolution of a system is described on a 
time loop (Keldysh contour, K) with "forward" ("— ") and 
'TDackward" ("+") directions (see Fig. [J), thus deflning 



two distinct branches for time evolution. Special care 
should be taken to deflne proper evolution operators for 
each direction. One could argue that the time evolution 
operator is the same along the full path, as in the Hermi- 

tian case, UK{t,to) = U+U^ = Tkc ^■^*o'^* (where 
Tk is the time-ordering operator along the Keldysh con- 
tour); which preserves the normaHzation along the con- 
tour, C/x(to ) = 1- This would mean that the time 
evolution would be governed by H on both branches. 
Even though this choice ensures that the "backward" 
evolution is the algebraic inverse of the "forward" one 
and it may seem a good choice, we should emphasize 
that it does not describe the physical process faith- 
fully. Physically, in the "backward" branch, the sys- 
tem should rewind its "forward" evolution and therefore 
should be the operator which dictates the dynamics. 
Then the "backward" time evolution operator is given 
by U+{to,t) = Ul{tM) = fe'^4'i*'^^(*') ^ Uz\t,t„), 
which is automatically anti-time-ordered (with T repre- 
senting the anti-time-ordering operator). Of course, one 
recovers the familiar result in the Hermitian case. We 
see that this approach would be hard to implement in 
a single-time formaHsm and a double-time approach is 
essential and rather natural for non-Hermitian problems. 
Thus, the "-I-" branch is governed by the Hermitian conju- 
gate theory and the time-ordered Green's function (G^~) 
is mapped into the anti-time-ordered one (G++) by Her- 
mitian conjugation; which is all natural in a time-loop 
formalism. This makes a generalized SK method an ap- 
pealing choice (30| . 
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FIG. 1: Sketch of the Keldysh contour. The arrows show the 
direction of the time evolution. Here G~~ is the time-ordered 
Green's function, G'*'^ is the anti-time-ordered Green's func- 
tion, is the "greater" Green's function and G""*" is the 
"lesser" Green's function. 
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To illustrate the formalism, we shall study a simple 
model of a single-mode asymmetric tunneling junction 
(see inset of Fig. [2]) j3l|, which captures essential aspects 
of irreversibility. We shall study the steady state trans- 
port of the system. The Hamiltonian of the model is 
given hy H — Hq + Ht, where 
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current operator in the usual way by / = ^dtAN, where 
AN = Nr — Nl, and Nl and are the particle number 
operators for the left and right leads, respectively. While 
the current is defined unambiguously for a symmetric 
junction (SJ), this is not the case for an asymmetric one. 
The above formula is incomplete if we do not specify on 
which branch AN resides. On the "— " branch, the cur- 



dx il}l{x)dxipr{x) (1) rent operator is given by / = \dtAN = ^[H,AN 



describes noninteracting left (L) a.nd right (i?) leads in 
a one-dimensional formulation [21|. iplix) {^r{x)\ rep- 
resent electron creation [annihilation] operators in the 
corresponding lead and v is the Fermi velocity which is 
taken to be the same for both leads. 

The irreversible nature of the problem is given by 
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which describes the tunneling with rates t*j^ in the 
asymmetric junction (A J) case. (In this letter, ks — h — 
e = l.) 

Let us derive a local theory which captures the essen- 
tial physics of the junction. After integrating out the 
leads exactly in a path-integral formulation and drop- 
ping the coordinate dependence of the electron creation 
and annihilation operators, the local action reads 



A 



'^{t)G-\t-t')'<i'{t')dtdt', 



(3) 



where ^' = ( V'l V'jj V'l V'fl ) • We have doubled the 
basis of the problem, which is customary in the SK for- 
mulation, and assigned extra indices "— " and to the 
fields according to the particular branch of the Keldysh 
contour they belong to. In Eq. we explicitly im- 
posed the steady-state constraint by specifying the time 
dependence of the inverse Green's function G~^; In the 
steady state, the Green's function depends only on the 
time difference, t — t' , and not on the "center-of-mass" 
time (i-t-f)/2 Thus, the Fourier transform (with 

respect to the time difference) is given by 
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where Sr = tanh '^2t'' ■ Here, we assumed that the leads 
are in equilibrium at their respective chemical potentials 
(/ii and /ifl) and temperatures (T^ and Tr). 

Using the Green's functions, we can calculate the 
transport properties of the asymmetric junction. 

Current - We first describe a current flowing across the 
junction. We expect that there should be an asymmetry 
between positive and negative voltage bias. We define the 



since H governs the time evolution on this branch. This 

prompts us to use the time-ordered (" " component) 

Green's functions for calculations. The current is thus 
given by 
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itLG],Z{cj)-tRGl]^{u;)]. (5) 



For consistency, we should get the same result if we 
do calculations on the branch. In this case, the cur- 
rent operator reads /+ = ^dtAN+ = AiV+]; in 
analogy with the "— " branch. Here W defines the time 
evolution and we should use the anti-time-ordered ("-I— 1-" 
component) Green's functions. In this case the current 
has the same form as in Eq. |[5]), but with time-ordered 
Green's functions replaced by anti-time-ordered ones and 
tL/R. ^ ^*R/L- Direct calculations indeed show that 
= = I. In the SJ limit, calculations can be done 
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FIG. 2: Current, /, as a function of voltage, V , for different 
temperatures (in arbitrary units). The current is shown for 
tn = tL = 0.05 and Ir = II ^ 0.05/^2 (SJ case, ioUed 
lines), for Ir = 0.05 and th ~ 0.05/\/2 (dashed lines) and 
for tn = 0.05 and tL — [solid lines). The temperature 
increases from bottom to top and Gq indicates the quantum 
of conductance. Inset: Sketch of the tunneling setup. 

analytically and the current is given by / = GV, where 
G ~ ^(]^^*^|2)2 f^S] and V — fiL — ^'R- The well known 
symmetry of the conductance in the SJ case, |t| 
is generalized to i^/L — ^ ^/^*l/r f^'" '^^^ 
results for the current in the weak tunneling limit are 
shown on Fig. [21 where we see the anticipated asymmet- 
ric behavior for unequal tunneling strengths. For large 
voltages the behavior is similar to the SJ case with the 
corresponding tunneling rate (t = tR,tL for positive or 
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negative bias, respectively). In the extreme limit when 
= 0, we see a clear diode-like behavior. The results 
also show temperature dependence (that disappears in 
the SJ Hmit), in particular, the current is nonzero for 
zero voltage due to finite temperature effects (/q oc T). 
For intermediate and large tunneling strengths, the cur- 
rent is significantly different from the SJ value even for 
large voltages. This can be seen on Fig. [31 which shows 
the conductance as a function of tunneling rate. 




FIG. 3: Large- voltage conductance as a function of tunneling 
strength (tjj) for a SJ, t = tR {dotted line) and in the extreme 
AJ case, tL = {solid line). 

The behavior of the conductance can be understood in 
terms of a perturbative expansion. For weak tunneling 
{t, ifl ^ 1), the expansions for the symmetric and asym- 
metric cases are identical {i.e., the leading orders coin- 
cide). While for intermediate tunneling (t, ifl ~ 1), many 
processes which are allowed for the S J are not allowed for 
the extreme asymmetric one, hence the reduction of the 
conductance. For large tunnehng strengths {t,tR ^ 1), 
the conductance decreases in the SJ case due to a reso- 
nance developing at the junction. In the AJ case, this 
resonance is suppressed and, therefore, the conductance 
increases, and saturates at the quantum of conductance 
in the limiting case (t^f — + oo). 

Noise - Now we turn to the noise, which is an in- 
dependent response [H, UHl- The current-current cor- 
relation function is usually defined in the following way 

S{t,t') = {:iit)::i{t'): + :iit')::i{t)-), (6) 

where : : denotes normal ordering. Bearing in mind the 
SK formaHsm, we redefine S as 

E {TK:i''{t)::r''{t'):), (7) 

which is equivalent to the definition in Eq. ([6]). 

From the definitions of the current operators on the 
corresponding branches, using Wick's theorem ji^ and 



Fourier transforming, we get the current noise power 

^ = tLtl{[G+ZoG],+]{iu) + [G],+ oG+Z]{cu)} 

+ tntU[G+], o GI+]{lo) + [GI+ o G+],]{lu)} 
~ tBtl{{Gl-i,oGll\{u,) + [GZA oG+^](n)} 
- ^Lt\{{G\l o G^+](c.) + [G^+ o G+^](c^)}, 

where the correlation product is defined as 

^G^{J ~uj)G-,{J). (8) 
In the SJ limit, we get the well established result [23 |: 

S'(tJ = 0) SttG^T + 2G(1 - 2txG)V coth ^ . (9) 

The dependence of ^(w = 0) on bias voltage and temper- 
ature, for weak tunnehng, is shown on Fig.lH The results 
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FIG. 4: Current noise power, S(aJ = 0), as a function of volt- 
age, V , for different temperatures (in arbitrary units). The 
line styles and parameters are the same as in Fig. [2] 

show similar asymmetric behavior as for the current. For 
large voltages, the behavior of the noise power is similar 
to that of the SJ with the corresponding tunneling rate 
(in exact analogy with the behavior of the current). In 
the Hmiting case with i_L = 0, we see the suppression of 
the noise for negative voltage bias, which is consistent 
with the diode-like behavior and is due to the reduction 
of the current in this regime (see Fig. [2]). In the shot noise 
limit {V > T), the Fano factor F = S{uo^Q)/{2I) = 1 
for week tunnehng as in the SJ case. 

Heat current - Now we shall turn to the thermal con- 
ductivity, calculate the heat current, Iq, and check the 
Wiedemann-Franz law. To define the heat current, we 
use the first law of thermodynamics to write 

(f AQ = dA£; + ^lLdNL - URdNR, (10) 

where AQ ^Qr-Ql and AE = (^Ho,r - Hq.l^. By us- 
ing the usual definitions for the currents, choosing a par- 
ticular branch ("— ") on the Keldysh contour and defin- 
ing the average chemical potential Ji, = {^il + M-r)/2, we 
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get Iq = Ie ~ P'I~ ■ Using the equations of motion for 
the electron creation and annihilation operators, and per- 
forming some algebra, the heat current reads 

/duj 
— {u;-mLGYL{^)-tRGl^{uj)]. (11) 

Our calculations on the "backward" branch yield Iq = 
Iq ~ Iq, in exact analogy with the result for the electric 
current. In the SJ case, we obtain the well known result 
m /q = if AT, with K = ^^+1$). , where AT = Tl - 
Tr and T = {Tr + Tl)/2. _The Wiedemann-Franz law 
holds and is given by Kj (GT) = 7r^/3 = Lq, the Lorentz 
number. 

In general, if bias voltages and temperature gradients 
are applied, the electric and heat currents are given by 

/ = LiiV + L12AT + lo, 
Iq = L21V + L22AT. (12) 

For the SJ case, Lu = G, L22 = K and L12 = £21 = 0; 
so no thermoelectric effects are observed. For the AJ 
case, expanding the result to linear order in V and AT, 
we see that still L12 — L21 — 0. (Notice that thermo- 
electric effects are observed in this case, but not at the 
linear order.) In terms of the transport coefficients, the 
Lorentz number is defined as L22/{LiiT). Our calcula- 
tions show deviations from the Wiedemann-Franz law. 
This behavior is due to the same effects as in the case of 
the conductance and the presence of the nonzero /o for 
the AJ case; Lq ~ tt^/S for small tunneling strengths, 
but it increases as the strength increases. 

We have presented a general scheme to deal with non- 
Hermitian Hamiltonians, and have shown that the gen- 
eralized Schwinger-Keldysh formalism is a natural choice 
for such applications. As an example, we have studied a 
simple model of an asymmetric tunneling junction where 
irreversibility is explicitly encoded into the Hamiltonian, 
hence making it non-Hermitian. We have calculated the 
steady state transport properties of the junction and seen 
that the behavior of the observables goes along with our 
physical expectations for the model. The approach is not 
limited to systems with steady states and can be used to 
address general non-equilibrium situations. We intend to 
apply it, for instance, to study the formation and evap- 
oration dynamics of interacting BECs, where a steady 
state is not achieved, and explore regimes {e.g., the uni- 
tary limit) that are not accessible with present methods. 
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